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Abstract 

The simultaneous determination of two unknown thermal coefficients for a semi-infinite 
material under a phase-change process with a mushy zone according to the Solomon- 
Wilson-Alexiades model is considered. The material is assumed to be initially liquid at 
its melting temperature and it is considered that the solidification process begins due 
to a heat ffux imposed at the fixed face. The associated free boundary value problem 
is overspecified with a convective boundary condition with the aim of the simultaneous 
determination of the temperature of the solid region, one of the two free boundaries of 
the mushy zone and two thermal coefficients among the latent heat by unit mass, the 
thermal conductivity, the mass density, the specific heat and the two coefficients that 
characterize the mushy zone. The another free boundary of the mushy zone, the bulk 
temperature and the heat flux and heat transfer coefficients at the fixed face are assumed 
to be known. According to the choice of the unknown thermal coefficients, fifteen phase- 
change problems arise. The study of all of them is presented and explicit formulae for 
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the unknowns are given, beside necessary and sufficient conditions on data in order to 
obtain them. Formulae for the unknown thermal coefficients, with their corresponding 
restrictions on data, are summarized in a table. 

Keywords: Phase Change, Convective Condition, Lame-Clapeyron-Stefan Problem, Mushy 
Zone, Solomon-Wilson-Alexiades Model, Unknown Thermal Coefficients. 

2010 AMS subjet classification: 35R35 - 35C06 - 80A22 

1 Introduction 

Heat transfer problems with a phase-change such as melting and freezing have been studied in 
the last century due to their wide scientific and technological applications [TlllllfJlIgl lT^lITTlfTHl 
[251135] . In particular, inverse problems related to the determination of thermal coefficients has 
been attracted many scientists because they are often ill-posed problems OTMMmZl 

In our recent work [7], we studied the determination of one unknown thermal coefficient for 
a semi-infinite material under a solidification process due to a heat flux imposed at the fixed 
boundary when the existence of a mushy zone it is considered following the Solomon-Wilson- 
Alexiades’ model [29] and the associated free boundary value problem is overspecified with a 
convective boundary condition [3T]. In this paper we study the simultaneous determination 
of two unknown thermal coefficients for the same phase-change process when some additional 
information on the physical phenomena is given. 

For the convenience of the reader, we repeat the main characteristics of the phase-change 
process considered in [3129]. The material is assumed to be initially liquid at a melting tem¬ 
perature of 0°C and the existence of three different regions in the solidification process is 
considered [29733] : 

1. liquid region at temperature T(x,t ) = 0: A = {(x,t) £ R 2 / x > r (t), t > 0}, 

2. solid region at temperature T(x,t ) <0: D s = {(x,i) G R 2 /0 < x < s(t), t > 0}, 

3. mushy region at temperature T(x,t) = 0: D p = {(x, t) G R 2 / s(t) < x < r(t), t > 0}, 
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where x = s(t ) and x = r(t) represent the free boundaries of the mushy zone, and T = T(x,t) 
represents the temperature of the material. The mushy zone is considered as isothermal and 
the following assumptions on its structure are made: 

1 . the material contains a fixed portion of the total latent heat per unit mass (see condition 
© below), 

2. its width is inversely proportional to the gradient of temperature (see condition ([5]) below). 

Finally, we recall that all of the thermal coefficients involved in the solidification process are 
assumed to be constants, where the bulk temperature — < 0 and the heat flux and heat 

transfer coefficients go > 0 and ho > 0 , respectively, at the fixed face are assumed to be known. 

In this paper, we consider that we also know the evolution in time of one of the two free 
boundaries of the mushy zone. More precisely, we assume that the free boundary x — s(t) has 
the form: 

s(t) = 2 cr y/t, t > 0 , ( 1 ) 

where a > 0 is a known coefficient. Thanks to this additional information on the physical 
phenomena, we will be able to determine simultaneously the temperature T = T(x,t ) of the 
solid region, the free boundary x = r(t) and two unknown thermal coefficients among the latent 
heat by unit mass l > 0 , the thermal conductivity k > 0 , the mass density p > 0 , the specific 
heat c > 0 , and the two coefficients 0 < e < 1 and 7 > 0 that characterize the mushy zone, by 
solving the following overspecified free boundary value problem: 


pcT t (x,t ) - kT xx (x,t) = 0 

0 < x < s(t), t > 0 

(2) 

T{s{t),t ) = 0 

t > 0 

(3) 

kT x (s(t),t) = pl[es(t ) + (1 - e)r(f)] 

t > 0 

(4) 

T x (s(t),t)(r(t ) - s(t)) = 7 

t > 0 

(5) 

r( 0 ) = 0 


( 6 ) 

kT x (0,t) = ^= 

t > 0 

(7) 

kT x (0,t) — h=(T(0, £) + DX) 

t > 0 

(8) 

Since inverse Stefan problems are usually ill posed, it is 

expected that restrictions on data 

has to be set in order to obtain solutions to problem ©-([HD- 

The goal of this paper 

is to obtain 
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necessary and sufficient conditions on data for the fifteen phase-change problems (l 2 i)- (l 8 |) that 
arise according to the choice of the unknown thermal coefficients, under which solutions can be 
obtained. Moreover, we also want to obtain those solutions explicitly. 

The organization of the paper is as follows: first (Sect. [2]) we prove a preliminary result 
where necessary and sufficient conditions on data for the phase-change process ([ 2 ])-® are given 
in order to obtain the temperature T = T(x, t) and the unknown free boundary x = r(t). Then 
(Sect. |3]), based on this preliminary result, we present and solve the fifteen different cases for 
the phase-change process ®)-®) corresponding to each possible choice of the two unknown 
thermal coefficients among l , k , p, c, e and 7 . Under certain restrictions on data, we prove that 
there are twelve cases where it is possible to fold a unique explicit solution and that there are 
infinite explicit solutions for the remainder three cases. At the end, explicit formulae for the 
unknown thermal coefficients, with their corresponding restrictions on data, are summarize in 
Table |T] 


2 Explicit solution to the phase-change process 

The following theorem represents the base on which the work of this article will be structured. 
Theorem 2.1. The solution to problem ®)-®) is given by: 


T(M) = -^e rtf’ 
k \ Ja 


erf 


2 \/od 


erf 


r(t) = 


y/c exp (cr 2 /a) 
Qo 


+ 2 a 


Vi 


0 < x < > 0 


t > 0 


(9) 

( 10 ) 


if and only if the physical parameters satisfy the following two equations: 


a + 


a 


y/c( 1 — e) exp (cr 2 /a) 


2<?o 


kD , 


Qo _ 

pi 

erf 

'a J q 0 y/air 

where the coefficient a, defined by: 

a = 


exp (cr 2 /a) 


00 1 ^ Qo 


h () Uqc: 


pc 


( 11 ) 

( 12 ) 

(13) 


is the thermal diffusivity. 
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Proof. The free boundary value problem (}2i)-(l8|) has the solution [29ll30lf33|| : 


T(x,t ) = A + B erf 



r(t) = 2 p\fod 


0 < x < s(t), t > 0 

(14) 

t > 0 

(15) 


where coefficients A, B and p have to be found. 

By imposing conditions ©-©, © and ([ 8 ]) on (Il4|i - f[l5|i we have that: 


A 




B 



and p 


7/0 exp (a 2 / a) 
2 q 0 \/a 


a 



(16) 


which corresponds to solution ©-(HDD, and that the physical parameters must satisfy equations 
(HU) and (J3U- ■ 


Hence, from Theorem 12.11 we have that there is an equivalence between solving the free 
boundary value problem ©-© with two unknown thermal coefficients and solving the system 
of equations (111H - 012D for the same two unknown thermal coefficients as in problem ©-©. 


3 Explicit formulae for the unknown thermal coefficients 


In this section we present and solve fifteen different cases for the phase-change process ©-© 
according to the choice of the two unknown thermal coefficients. With the aim of organizing 
our work, we classify each problem by making reference to the coefficients which is necessary 
to know in order to solve it (see Theorem ED: 

Case 1: Determination of e and 7 , Case 2: Determination of e and l , 

Case 3: Determination of 7 and l, Case 4: Determination of e and k , 


Case 5: Determination of e and p, 
Case 7: Determination of 7 and k , 
Case 9: Determination of 7 and c, 


Case 6 : Determination of e and c, 
Case 8 : Determination 7 and p, 
Case 10: Determination of l and k , 


Case 11: Determination of l and p, Case 12: Determination of l and c, 
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Case 13: Determination of k and p, Case 14: Determination of k and c, 
and 


Case 15: Determination of p and c. 


Moreover, we introduce several functions and parameters which are named with an index 
according to the number of the case where they arise for the first time. 


Theorem 3.1 (Case 1: determination of e and 7 ). If we consider the phase-change process (|2jl- 
() 8 P with unknown thermal coefficients e and 7 , then it has infinite solutions given by m-m 
with: 

7 = k(i°-e) {~pt 6XP “ X ) exp ( -a2 /a) ( 17 ) 

and any e G (0,1), if and only if the remainder physical parameters satisfy condition (fT2l) and 
the following inequality: 

0 < — 7 — exp (—cr 2 / ol) — 1. (Rl) 

pier 

Proof. Due to Theorem l2.ll we have that the phase-change process (|2|) -(l8|) has the solution given 
in (l9D- (fT0D if and only if e and 7 satisfy equation (fill) and the remainder physical parameters 
satisfy condition f| 12 |) . Then, we have from equation (fill) that 7 must be given by (1T7D for any 
e G (0,1). To finish the proof, only remains to observe that this coefficient 7 is positive if and 
only if inequality (IRlj) holds. ■ 


Theorem 3.2 (Case 2: determination of e and l). If we consider the phase-change process (121) - 
P) with unknown thermal coefficients e and l, then it has infinite solutions given by ©-(HDD 
with: 


1 = 


g 0 exp (—a 2 / a) 


(18) 


P° 1 + ex P ( a 7«) 

and any e G (0,1), if and only if the physical parameters h 0 , q 0 , D^, a, p, c and k satisfy 
condition (Tl 2 l) . 


Proof. It is similar to the proof of Theorem 13.11 


Theorem 3.3 (Case 3: determination of 7 and l). If we consider the phase-change process 
©-© with unknown thermal coefficients 7 and /, then it has infinite solutions given by ©- (1IU1) 
with any 7 > 0 and l given by (TT 8 lh if and only if the parameters h 0 , q 0 , D 0c , a, p, c and k 
satisfy condition (Tl 2 l) . 
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Proof. It is similar to the proof of Theorem 13.11 


Remark 1. Let us observe that it follows from the previous three theorems that, under certain 
conditions for the data of the problem, the phase-change process (J2])-(]S]) corresponding to Cases 
1, 2 and 3 has an infinite number of solutions. 

Theorem 3.4 (Case 4: determination of e and k ). If we consider the phase-change process 
P|)- (18j) with unknown thermal coefficients e and k, then it has the solution given by (lUjb flTH 
with: 


6 = i - m) 


k = pc 


where f is the only positive solution to the equation: 

apcD oo ( q 0 
g 4 {x) = -— I 1 - 


QoV^ V h o D oo 

and the real functions / 4 and g 4 are defined by: 


/ 4 (x) = —— ( -^-exp (— x 2 ) — 1 ) x 2 exp (— x 2 ) and g^x) = a;erf (x), x > 0, 


7 pc<T \pla 

if and only if the remainder physical parameters satisfy the next three inequalities: 


0 < 1 - 
0 < 1 - 
<?o 


Qo 


hoDoo 

%_ 

pla 


1 - 


< -- 9 a 


hoDoa apcD. 

and any of the following three groups of conditions: 
Group 1: 


'In 


Qo_ 

pla 


/ 4 (h) > 1 

9o 




or 


<?o 


> -“—34(C2) 


/ioDco opcDm hnDtt ' apcD c _ 

where ()i and (2 are the only two positive solutions to the equation: 


(19) 

( 20 ) 


(E4) 


( 21 ) 


(R2) 

(R3) 

(R4) 


(R5) 

(R6) 


f,(x) = 1, 


( 22 ) 
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and t] is the only positive solution to the equation: 


Group 2: 


% 

pier 


(1 - 2x 2 ) 


(1 — x 2 ) exp(x 2 ). 


Group 3: 


h(v) = 1 

-I Qo , qoy/n , \ 

1 _ nr ^ 

IiqJ^oo crpcDoo 

where p is the only positive solution to equation (1231) . 

Mv) < 1 


where 7 ] is the only positive solution to equation (l23|) . 


( 23 ) 


(R7) 

(R8) 


(R9) 


Proof. As in the proof of the Theorem 13.11 we have from Theorem 12.11 that the phase-change 
process o-« has the solution given in (l9|i- iTTCTI) if and only if e and k satisfy equations (TIT]) 
and (fT2j) . Introducing the following dimensionless parameter: 


£ = 



(24) 


we have that the solution of the system of equations (TTT1) - (1T2T) is given by (fimhfliZUj) if and only 
if f is a solution to equation (IE4I) . Then, we need to prove that the restrictions on data given 
in the statement are necessary and sufficient conditions for the existence of a positive solution 
to equation (1E4|) and for obtain that the coefficient e given in (Il9]i is a number between 0 and 
1. 

We first note that equation (1E4I) admits a positive solution if and only if inequality (1R 21) 
holds, because g 4 is an increasing function from 0 to +00 in M + . Henceforth, we will assume 
that (1R2D holds. 

Let us now focus on the fact that e G (0,1). On one side, we have that e is less than 1 if 
and only if (see 03): 

which is equivalent to inequality (IR3I) and: 


f < log 



(25) 






Since c/ 4 is an increasing function and £ satisfies equation (IE4D . applying function g 4 side by 
side of inequality (1251) we have that it is equivalent to inequality (IR4I) . Therefore, from now on 
we will assume that inequalities (1R 3[) and (1R4]) also hold. 

On the other side, we have that e is positive if and only if (see 03): 

MS) < 1. (26) 


Let us now think of the existence of extreme values of / 4 . We have that / 4 (x) = 0 if and only if 
x is a positive solution to equation (1251) . In order to study equation (125]) . let us introduce the 
real functions w 4 and u 4 defined by: 


% 


u*(x) = ^-( l- 2 x 2 ) 
pier 


and v±(x) = (1 — x 2 ) exp (x 2 ), x > 0. 
Since w 4 and n 4 are decreasing functions such that: 


(27) 


w 4 (0 + ) = > 1, w 4 (l/\/2) = 0, m 4 (+oo) = 

pi a 

u 4 (0 + ) = 1, n 4 (l) = 0, 


-oo, 


u 4 (+oo) = —oo, 


it follows that equation (125]) admits only one positive solution p. Hence, / 4 admits only one 
critical point, p. Moreover, / 4 (x) > 0 if and only if 0 < x < p. Then, / 4 has a relative maximum 
in p. Since / 4 (0 + ) = 0 and / 4 (+oo) = 0, it follows that / 4 reaches its absolute maximum in p 
and that / 4 (r/) > 0. 

Now, we will study three different situations: f(rj) > 1, f(rj) = 1 and f(p) < 1, which are 
related to conditions given in Group 1, Group 2 and Group 3, respectively. 

If f(rj) > 1, that is, if inequality (IR5I) holds, we have that £ satisfies inequality (125]) if and only 
if: 


f < Ci or £ > £ 2 , 


(28) 


where £ 4 and £ 2 are the only two positive solutions to equation (122]) . Applying the increasing 
function g 4 side by side to both inequalities and taking into account that £ satisfies equation 
(]E4[) . it follows that (128]) is equivalent to (1R 61) . 

If f( 7 l) = I; that is, if (IR7I) holds, we have that £ satisfies inequality (125]) if and only if £ 7^ p. 
We now proceed as in the previous situation and obtain that £ 7^ p is equivalent to (IR8j) . 
Finally, if f(p) < 1, that is, if (1R9D holds, we have that inequality (125]) holds immediately. ■ 
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The previous Theorem 13.41 states necessary and sufficient conditions for the data in the over- 
specified free boundary value problem ®-® under which it is possible to fold the temperature 
T = T(x,t), the free boundary x — r(t) and the two unknown thermal coefficients e and k. 
Nevertheless, there are also some sufficient conditions on data, which are easier to check than 
the necessary and sufficient conditions given in Theorem 13.31 that enable us to fold the solution 
to problem ([2])-®. Next Proposition is related to those sufficient conditions. 


Proposition 3.1 (Sufficient conditions for Case 4). Let problem with unknown ther¬ 

mal coefficients e and k. If the remainder physical parameters satisfy inequality Hi and the 
following three conditions: 


qoVn 

apcD c 


-9 4 


0 < 


2gp 

P'yca 



< 1 - 


go 

h() Coc 


< 


qoVn 

-n -5,4 

apcD^ 




(RIO) 

(RH) 


where: 


then the solution to problem (E))-([I|i is given by m-m with e and k given by m an d ( PSj) ; 
being f the only positive solution to the equation m- 



pl<r , ( 9o 

v 4 = -In —— 

2 g 0 \plcr 


1 + 


1 + 


2 7c 


\ '>"(&) 


Proof. Let us assume that inequalities (IR3I) , OR 101) and 0R11I) hold. We have seen in the proof 
of Theorem 13.41 that e and k must be given by (1T9D and (12UD . and that f must satisfy equation 
0E4j) . We have also seen that (IR2|) is a necessary and sufficient condition for the existence and 
uniqueness of the solution to equation (1E4|) . Then, since the first inequality in 0R1Q|) implies 
OR 21) . we have that equation 0E4I) admits only one positive solution. Moreover, we have seen 
that e given in (lT9l) is less than 1 if and only if inequalities 0R3j) and (1R4|) hold. Since inequality 
(1R10I) implies inequalities (1R 31) and (IR4j) . we have that e given in (TT9]) is less than 1. Finally, 
we have also seen that e given in (1191) is positive if and only if inequality (129|) holds. The rest 
of the proof will be devoted to demonstrate that the first inequality in (IR10D and (IR11I) imply 
inequality (126]) . 

Since e given in (fT9j) is positive, we have that ^ exp (—£ 2 ) — 1 > 0, that is: 
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Then we have (see (l2lil ): 

(30) 

From the above, it follows that it enough to prove that the first inequality in (IR10D and (IRllj) 
imply that the right hand side of (l30l) is less than 1. Let wq the function defined by: 


uq(x) = a 4 x 2 — b^x — 1, x > 0, 


( 31 ) 


with aq and 64 given by: 

2 q 2 o 


a 4 = 


In 


<?o 


> 0 and 64 = —— In () > 0. 

pajc \pzoq 


(32) 


p 2 o 2 l^c \pl<7 , 

We have that z/ 4 given in (129|) is a positive root of uq. Moreover, we have that inequalities (1R3I) 
and (IR11D imply u 4 < 1. Since the another root of «q is negative, we have that the right hand 
side of (130| is less than 1 if and only if exp (—£ 2 ) < z/4, that is, if and only if: 


5 > ,/ln I - 


Only remains to observe that this last inequality is equivalent to the first inequality in (IRlOp 
because is an increasing function and £ satisfies equation 


Theorem 3.5 (Case 5: determination of e and p). If we consider the phase-change process 
Q-dSD with unknown thermal coefficients e and p, then it has the solution given by (l91) - (fTUl) 
with: 


e= Wb(0 

k fC 2 

P = -c{a 

where £ is the only positive solution to the equation: 


9s(x) = 


kDoo f 


q 0 cry/7T V h 0 D o 
and the real functions /s and are defined by: 


Qo 


(33) 

(34) 

(E5) 


erf (x) 


/5(a) = ^7 eXP ^2 X ~ _ ex P ( -a;2 ) and 95(2) = —7-^) X > 0, ( 35 ) 


if and only if the remainder physical parameters satisfy the following two conditions: 


gocq/i 

kD c 


35 (Ci) < 1 


Qo 


ho D 


o-r/00 


. f 2 q 0 a qocr^TT , 

<mm iHtI’W 95(C2 ' 


(R12) 
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where Ci and ( 2 are, respectively, the only positive solutions to equations: 


qoac 

Ik 


exp (— x 2 ) 


= x 


2 


q 0 ac 

Ik 


exp (— x 2 ) 


7 k 
2q 0 a 


exp ( x 2 ) + 1 


x 


(36) 

(37) 


Proof. As in the previous proofs, we have from Theorem 12.11 that the phase-change process 
(12i)-Qgj) has the solution given in (15jh (lTUD if and only if e and p are given by (1551) and (J2H), 
where the dimensionless parameter f (see (OJ)) is a solution to equation (IE 51) . Then, we need 
to prove that the two restrictions given by (1R12|) are necessary and sufficient conditions for the 
existence of a positive solution to equation (IE5j) and for obtain that the coefficient e given in 
(mi) is a number between 0 and 1. 

Since g$ is a decreasing function from to 0 in R + , we have that equation (1E5|) admits 
positive solutions if and only if: 


0 < 1 - 


go 


2g 0 a 

kDoo 


(38) 


Let us assume for a moment that (138|) holds and focus on the fact that e e (0,1). As we 
did in the proof of the Theorem 13.41 let us first think about the existence of extreme values of 
f 5 . We have that f$(x) = 0 if and only if: 


u 5 (x) = x 2 , 


(39) 


where is the real function defined by: 


u 5 (x) = 


q 0 crc 

Ik 


x* 


1 + — ) exp ( -x 2 ), X > 0. 


(40) 


Since u§ is an increasing function from 0 to +oo in M + , it follows that equation (l39j) admits 
only one positive solution q. Then f 5 has only one critical point, q. Moreover, /^(x) > 0 if and 
only if x > q. Then, / 5 has a relative minimum in q. Since /s(0 + ) = +oo and /5(+oo) = 0, 
it follows that /s reaches its absolute minimum in q and that /s(^) < 0. Then we have that 
the coefficient e given in (1331) is a number between 0 and 1 if and only if £1 < £ < C2, where <Ci 
and (2 are the only two positive numbers such that /s(Ci) = 0 and / 5 (( 2 ) = 1, that is, the only 
positive solutions to equations (IMll and (1571) , respectively. Since g 4 is a decreasing function and 
£ satisfies equation (lE5jh we have that Ci < £ < C2 is equivalent to: 


qoOjs/TT 

kD r 


95((i) < 1 


qo 


hc)D, 


O-^oo 


^ qo°Vn ,,, 

K hvT 95(Ci) ' 


(41) 
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Only remains to observe that inequalities (138D and (HTT) are equivalent to the inequalities given 
by (EED. ■ 

Theorem 3.6 (Case 6 : determination of e and c). If we consider the phase-change process 
o-o with unknown thermal coefficients e and c, then it has the solution given by m-m 
with: 


£ = i - /»(£) 



(42) 

(43) 


where £ is the only positive solution to the equation flE5jl . and f e is the real function defined 


by: 


fe(x) = ~y (jy exp ) exp ( _a:2 )’ x > °> 


if and only if the remainder physical parameters satisfy the following condition: 

q 0 a^/7T 


kD r 


-h Win 


%_ 

pla 


< 1 - 


% 

HqDqq 


and any of the following two groups of conditions: 
Group 1: 


-*L>Uk + 1 

pla 2q 0 a 


1 - 


% 


hoDoo 
where: 


< mm 


2 goer 

k Dcxd kD cxd 


A I'/>■>(- 


^6 = 


pla 

2 g 0 


i + W i + 


2'yk 

a 2 pl 


Group 2 : 




g 0 2 q 0 a 


kD qq 

Proof. It is similar to the proof of Theorem 13.41 


(44) 

(R13) 

(R14) 

(R15) 

(45) 

(R16) 

(R17) 


Theorem 3.7 (Case 7: determination of 7 and k). If we consider the phase-change process 
(ED-dHD with unknown thermal coefficients 7 and k, then it has the solution given by ©-(HDD 
with: 


7 = 


2 go 

apc(l — e) 



exp (—£ 2 ) 



£ 2 exp (—£ 2 ) 


(46) 
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and k given by 020]) . where £ is the only positive solution to the equation 0E4I) . if and only if 
the physical parameters h 0 , g 0 , D 00 , cr, /, p and c satisfy conditions 0R2D . 0R3I) and OR 41) . 


Proof. Once again, we have from the Theorem 12.11 that 7 and k must be given by (1451) and (1201) . 
where the dimensionless parameter £ (see (1241) 1 is a solution to equation (1E4D . As we saw in the 
proof of the Theorem 13.41 the equation (1E4|) admits positive solutions if and only if inequality 
(1R.2D holds. 

To complete the proof only remains to observe that the coefficient 7 given in (1401) is positive 
if and only if 0 < ^exp(—£ 2 ) — 1 which, as we also saw in the proof of Theorem 13.41 is 
equivalent to inequalities (1R3() and (1R 41) . ■ 


Theorem 3.8 (Case 8 : determination of 7 and p ). If we consider the phase-change process 
(12])- (l8|) with unknown thermal coefficients 7 and p, then it has the solution given by (l9])- (fT0l) 
with: 


7 = 


(47) 


2 qua (®or exp (-£ 2 ) \ , 

k(l-e)\lk £ Ijexpm 

and p given by (]34|) . where £ is the only positive solution to the equation (1E5D . if and only if 
the physical parameters h 0 , g 0 , D 0c , a, k and c satisfy conditions OR 21) . OR 171) and: 


kD oo / g 0 \ 

95 9 qocry/n \ hoDoo) 

where g 5 is defined in 040]) and 7 is the only positive solution to the equation: 

goccx exp(—x 2 ) 

Ik x 


(R18) 


(48) 


Proof. From the Theorem 12.11 we have that 7 and p must be given by 0471) and (154|) . where 
the dimensionless parameter £ (see (124]) ) is a solution to equation (IE5j) . As we saw in the proof 
of the Theorem 13.51 the equation 0E5D admits a positive solution if and only if inequality (138|) 
holds, that is, if and only if inequalities (IR2[) and OR 171) hold. 

We also have that the coefficient 7 given in (1471) is positive if and only if: 


m 8 (£) > 0, 


(49) 


where u$ is the real function defined by: 


u 8 (x) 


Qoco exp (-x 2 ) _ 
Ik x 


x > 0 . 


(50) 
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Since % is a decreasing function from +oo to —1 in M + , we have that equation (l48]i has only 
one positive solution p. Therefore, inequality (E© holds if and only if £ < rj. Only remains to 
observe that £ < rj is equivalent to condition (1R.18|) because r/ 5 is a decreasing function and £ 
satisfies equation (1E5D . ■ 

Theorem 3.9 (Case 9: determination of 7 and c). If we consider the phase-change process 
©-« with unknown thermal coefficients 7 and c, then it has the solution given by m-m 
with: 


7 = 


2 q 0 cr 


Qo 


exp (-a -1 exp (-a 


(51) 


k( 1 — e) 

and c given by (l43lh where £ is the only positive solution to the equation (1E5|) . if and only if 
the physical parameters h 0 , q 0 , a, l, k and p satisfy conditions (IR3D . (IR13|) and (IR17D . 


Proof. It is similar to the proof of Theorem 13.81 ■ 

Theorem 3.10 (Case 10: determination of l and k ). If we consider the phase-change process 
©-« with unknown thermal coefficients l and k, then it has the solution given by ©-(HDD 
with: 


l = 


Qo_ 

pa 


pca^(l—e) 


1 + 


2<jo 


exp (£ 2 

£ 2 


exp (—£- 


(52) 


and k given by (l2T)]h where £ is the only positive solution to the equation (1E4H . if and only if 
the physical parameters h 0 , q 0 and D 00 satisfy condition (1R2[) . 

Proof. It is similar to the proof of Theorem 13.71 ■ 

Theorem 3.11 (Case 11: determination of l and p). If we consider the phase-change process 
(ED-dHD with unknown thermal coefficients l and p, then it has the solution given by (l9]h (fT0]i 
with: 


l = 


q 0 ca 

k 


1 + 


7&(1—e) 


exp (£ 2 


exp(—£ 2 ) 

£ 


(53) 


and p given by (l34]h where £ is the only positive solution to the equation (lE5|h if and only if 
the physical parameters h 0 , q 0 , D^ and k satisfy conditions flR2|) and (IR17[) . 


Proof. It is similar to the proof of Theorem 13.81 


Theorem 3.12 (Case 12: determination of l and c). If we consider the phase-change process 
m-m with unknown thermal coefficients l and c, then it has the solution given by ®-(Ho]) 
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with: 


i= ®L 

pa 


1 + ^exp (^) 


exp (-£ 2 ) 


(54) 


and c given by (1431) . where £ is the only positive solution to the equation (1E5j) . if and only if 
the physical parameters h 0 , q 0 , and k satisfy conditions (1R 21) and (IR17D . 

Proof. It is similar to the proof of Theorem 13.81 ■ 

Theorem 3.13 (Case 13: determination of k and p ). If we consider the phase-change process 
®-o with unknown thermal coefficients k and p, then it has the solution given by (]9]) - (fl0]) 
with: 

q 0 ay/TT 


k = 


n [1_2s_ 

1 1 h 0 D o 


&>(£) 


P = :F^- 04 (f) 


(55) 


(56) 


caDoo 

where real functions ,ry 4 and g 5 are defined in (I2TT) and (1551) . £ is the only positive solution to 
the equation: 

9 ia{x) = —hi 3 (x), (E13) 

C13 

and the real functions (713 and h\ 3 are defined by: 


exp (—x 2 ) 

gi 3 (x) = - — —r— and hi 3 (x) = x + 613 exp ( x 2 ) erf (x), x > 0, (57) 


erf ( x ) 


with: 


2cRoo / q 0 


ly/n 


O-^oo 


6i , = 7^(1-«) <?» 


9 D (1_ q2 ~ 

1 1 h 0 D o 


h n D, 


O-^oo 


(58) 


if and only if the physical parameters ho, g 0 and D^ satisfy condition 


Proof. We have from Theorem 12.II that the phase-change process ([5])-(f8]) has the solution given 
in (EJD-dTUD if and only if k and p satisfy equations dll]) and (115]) . Since inequality OR 21) is a 
necessary condition for the existence of a solution to equation (115]) . now and on we will assume 
that inequality (1R2j) holds. 

We have that the solution to the system of equations (TTTI) - (TT5]) is given by (1551) and (156]) . 
where the dimensionless parameter £ is a solution to equation (]E13j) . Only remains to observe 
that equation 0E13j) admits a positive solution since g\ 3 is a decreasing function from +00 to 0 
in M + and h \3 is an increasing function from 0 to +00 in M + . ■ 
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Theorem 3.14 (Case 14: determination of k and c). If we consider the phase-change process 
((2))-© with unknown thermal coefficients k and c, then it has the solution given by ©- (fill 
with k given by (1551) and: 


c = 


qoV^ 


apD , 


1 _ go 

00 ' - 1 - h 0 D o 


94®, 


(59) 


where the real function c/4 is defined in f|2T]) . £ is the only positive solution to the equation: 

ai40u (®) = h u (x) (E14) 

and the real functions c/14 and hu are defined by: 

gu(x) = ( —p- exp (— x 2 ) — 1 ) x and h\Ax) = erf (x) exp (x 2 ), x > 0, (60) 

\pla J 

with: 


2 D r 


a 14 — 


0 o 


7\/7r(l -e) h 0 D OO/ 
if and only if the remainder physical parameters satisfy conditions (1R2|) , (1R3() and: 

9u(v) > h 14 (ri), 

where r) is the only positive solution to the equation: 


Qo_ 

pla 


(1 - 2x 2 ) = exp ( x 2 ). 


(61) 


(R19) 


(62) 


Proof. It is similar to the proof of Theorem 13.131 ■ 

Theorem 3.15 (Case 15: determination of p and c). If we consider the phase-change process 
(J2 |) - (l8l) with unknown thermal coefficients p and c, then it has the solution given by ( © -([TO ]) 
with: 

q 0 exp (—£ 2 ) 


P = 


la 1 + ex P (^ 2 ) 

kl 


c = 


vqo 


1 + 


2q 0 a 

7/0(1 — e) 


2 q 0 cr 


exp (f 2 


^exp (£ 2 


(63) 


(64) 


where £ is the only positive solution to the equation (1E5|) . if and only if the physical parameters 
h 0 , q 0 , Doo, a and k satisfy conditions (IR2I) and OR 171) . 

Proof. It is similar to the proof of Theorem 13.131 ■ 

Next table summarizes the results of this section corresponding to the fifteen cases for the 
simultaneous determination of two unknown thermal coefficients for problem (|2|) -(]8|). 
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Case Thermal coefficients Equation for Restrictions on data 



1 0 < e < 1 

7 = ]|SSj (& exp (-aya) - l) exp (-<x 2 /a) 

2 0 < e < 1 

J _ go exp (—o- 2 /g) 

P fJ [ 1 + 7 2 ( g 1 0 7 ) exp (o~ 2 /a)] 

3 7 > 0 

J _ gocxp(-o- 2 /a) 

[ :1 + 2 ^T' exp (o - 2 /a) J 

4 e = 1 - / 4 (0, with U given in (ED B 

k = pc (f ) 


5 e = 1 -/ 5 (0, with/ 5 given in (J25D dE5D 

p = !(f) 2 

6 e = 1 - / 6 (0, with f 6 given in (01 (El 


^ 7 = ($ exp (-e) - l) f 2 exp (-f 2 ) » 

k = pe (|j 

8 7 = ^5 (T ^ - 1 ) exp (-J 2 ) (E 3 

P = !(f) 2 


0 , CD 


ca 


ca 


a, ca, ca, ca, ca 

or 

dm , dm , dm , dm, m 

or 

dm, dm, dm, dm 

dEH 


dm, dEUD, dra 

or 

dm , (ThTITTi) , dEUD 

dm- dm- dm 


dm , dnizD- dM 
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Case Thermal coefficients Equation for Restrictions on data 



7 = ^y feex pi 


- 1 exp (-f 


p \ a ) 


, dEUD, dE32D 


10 


l = 33. 

pa 


1 exp (g 2 ) 

l | peg7(1 —e) £ 2 


k = pc (f) 


2 10 


exp (-£ 2 


m 


dEI 


ii 


7 _ go ca 
k 




exp (-g 2 ) 
( 


p = HiY 


12 


l = 33. 

pa 


!+^a“P(C) 


exp (-£ 2 ) 


p \ a / 


13 k = 


qoay/Tr 


Doo 1- 


h 0 D o 


95 (0> with g 5 given in 


P = f^74(0> with ^4 given in (ED 


(EH 


», dEUD 


dE2D, dEUD 


dE2D 


14 k= n /T^o 1 ^(0) with S iven in (ED 
°° l h 0 D °o7 

-^4(0) with 5(4 given in (ED 


c = 


gov^ 




dm 


dE2D, dEHD, dEED 


15 


p = 20 


exp(—g 2 ) 


^ i+^a- PK 2 ) 


C = 


fcZ 

o-go 


f 2 exp(f 2 ) 


dEED 


dE2D, dEED 
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Case 

Thermal coefficients 

Equation for 

£ = = ffTf 

^ yj OL V k 

Restrictions on data 



Table 1: FORMULAE FOR PROBLEM (121)- (!51) . Explicit formulae for the two unknown thermal 
coefficients among l , 7 , e, k , p or c, with the corresponding equation for the dimensionless 
parameter £ (see flMll ) and restrictions on data. 


4 Conclusions 

In this paper, we consider a semi-infinite material under a solidification process caused by 
an initial heat flux boundary condition, when the thermophysical parameters involved in the 
phase-change process are assumed to be constant and the Solomon-Wilson-Alexiades’ mushy 
zone model it is considered. We suppose that evolution in time of one of the two free boundaries 
of the mushy zone, the bulk temperature and the heat flux and heat transfer coefficients at 
the fixed face are known. We overspecify the associated free boundary value problem with 
the aim of the simultaneous determination of the temperature of the material, the unknown 
free boundary and two unknown thermal coefficients among the latent heat by unit mass, 
the thermal conductivity, the mass density, the specific heat and the two coefficients that 
characterize the mushy zone. We first prove a preliminary result where necessary and sufficient 
conditions on data for the phase-change process are given in order to obtain the temperature 
and the unknown free boundary. Then, based on this preliminary result, we present and solve 
the fifteen different cases for the phase-change process corresponding to each possible choice 
of the two unknown thermal coefficients. We prove that, under certain restrictions on data, 
there are twelve cases where it is possible to find a unique explicit solution and that there are 
infinite explicit solutions for the remainder three cases. For each case, we give formulae for the 
temperature, the unknown free boundary and the two unknown thermal coefficients, beside the 
necessary and sufficient conditions on data in order to obtain them. We summarize explicit 
formulae for the two unknown thermal coefficients in Table [E] 
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Appendix: Summary of restrictions on data, definition of 
functions and equations for £. 

List of restrictions on data 


0 < exp (— a 2 /a) — 1 

pl(J 


0 < 1 


Qo 


hnD 


O-'-'oo 


0<^-l 


pla 
% 


< - —9i 


HqDoo ' apcD 

h(v) > 1 



where rj is the only positive solution to the equation: 
-y-(l — 2x 2 ) — (1 — x 2 ) exp(x 2 ), x > 0 


(Rl) 

(R2) 

(R3) 

(R4) 

(R5) 


«o 


<-^-S4(Ci! 


or 


Qo 


^ Qoy/Tt / /- \ 
>-- 9a\Q2) 


hoD^ apcDoo hoD^ ' apcD c 

where £i and ( 2 are the only two positive solutions to the equation: 


(R6) 


f 4 (x) = 1, x > 0 


Mv) = 1 


where rj is the only positive solution to the equation (123|) . 
Qo , QoV^ 




-9«kn) 


hoDoo ' apcD c 

where rj is the only positive solution to the equation 

Mv) < 1 

where rj is the only positive solution to the equation (123|) 
qoy/n 


apcD c 


-9a 


1 In ( — 


< 1 - 


go 


^ qoV^ 

< - T^9a 

(J pcDoo 


'In 


go 

pier 


(R7) 


(R8) 


(R9) 


(RIO) 
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where: 


^ In ( 

2 go VpW 


1 + 



27c 

,ln (^) 


2g 0 


0 < 

P'yca 

qo& 


In 


_qo_\ f go 

pi a J \pl <7 


kD r 


#5(Ci) < 1 — 


go 


hoD r _ 


< mm 


2g 0 q~ %ojA 
’ kD r 


05 (C2 


(RH) 

(R12) 


where Ci and C2 are, respectively, the unique only solutions to equations: 


%ac 

Ik 

%ac 

Ik 


exp (—x 2 ) = x 2 , x > 0 


exp (—x ) = 


7 k 
2 g 0 a 


exp (x 2 ) + 1 


x 


x > 0 


-9 5 
7/c 


goo^/n 

kDoo 

— > _ 

p/tr 2<7 0 ct 

hoDoo 
where 


'In [ 


plo 


< 1 


Zo 


h(\D, 


O-^oo 


+ 1 


< nun 


Zqocr gpcry/n 

kD^ ’ kD 

pi a 


-95 


'In 


^6 = 


2<?o 


1 + 1 /1 + 


27 k 

cr 2 pi 




<?o 


hnDr, 


< 


2 g 0 a 
kD r 


L 0 - Ly oo ,XjJ - y oc> 

M <22^ f 1 - 


Zo 


goO’v^ 7 V h 0 D c 

where 77 is the only positive solution to the equation: 

goca exp(—x 2 ^ 

Ik x 


= 1, x > 0 


(R 13 ) 

(R 14 ) 

(R 15 ) 


(R 16 ) 

(R 17 ) 

(R 18 ) 


5-14(77) > hu(r]) (R 19 ) 

where 77 is the only positive solution to the equation: 

-^-(1 — 2x 2 ) = exp (x 2 ), x > 0 
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List of definitions of functions (x > 0) 


2 g 0 ( go 


fa(x) = 
fa(x) = 

9is(x) = 


- ! 

7 pea 

iyp/a 

2 q 0 a i 

' qoccr 

7 /c \ 

v Ik 

2g 0 (T / 

' go c 

7 A: \ 

v p/cr 

exp (- 

-a; 2 ) 


-1 exp (-a; 2 ), g 5 (x) = 


erf (x) 


x 


erf (x) 


hi3(x) = x + 613 exp ( x 2 ) erf (x), 


with: 


( q 0 
« 13 = — (i - — 


IsJk 


O-^oo 


; 7 V^(l - e) 0 / <?o 

Ol3 — -7-T) c 13 — 21 — 




hnD 


O-^oo 


0i4(z) = ( exp ( -X 2 ) - 1 ) x, 


/iu(x) = erf (x) exp (x 2 ), 


with: 


2 D r 


a i4 — 


7V^(1 - e ) V h 0 D, 


1 - 


Qo 


List of equations for £ 


34 (x) 


Q-pcPpo f _ go 

qoV^ V hoDoo 


95{x) : 

913{x) 


kD r 
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g o o-0r V 

a 13, / s 

-«13W 

Cl3 


go \ 
hoDoo ) 


au9u(x) = h u (x) 


(E4) 

(E5) 

(E13) 

(E14) 
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